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Calculations of the breakup of 8 B and 11 Be are performed with the aim of analyzing their sensi- 
tivity to the projectile description. Several potentials adjusted on the same experimental data are 
used for each projectile. The results vary significantly with the potential choice, and this sensitivity 
differs from one projectile to the other. In the 8 B case, the breakup cross section is approximately 
scaled by the asymptotic normalization coefficient of the initial bound state (ANC). For n Be, the 
overall normalization of the breakup cross section is no longer solely determined by the ANC. The 
partial waves describing the continuum are found to play a significant role in this variation, as the 
sensitivity of the phase shifts to the projectile description changes with the physical constraints 
' imposed to the potential. 
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I. INTRODUCTION 

Since the early days of radioactive beam experiments, breakup reactions have been an important source of infor- 
mation on the structure of nuclei near the dripline. Q, 0, S 0, ■ Usually, the exotic nucleus — the projectile — is 
simulated by a loosely bound two-body system. The spectroscopic factor of that single-particle wave function is 
obtained from the comparison of the theoretical predictions with the data [S • Alternatively, it has been 
suggested 0, that asymptotic normalization coefficients (ANC) can be obtained from breakup measurements. 
\Q , The basic idea therein is that breakup reactions of loosely bound nuclei are highly peripheral. Since the asymptotic 
behavior of the projectile wave function is in general well known, the breakup cross section is proportional to the 
square of the ANC. In this work, we examine the validity of such procedures. We focus on the dependence of the 
■ normalization of the Coulomb induced breakup cross section on the single particle potential used for the description 
of the ground state, as well as the dependence on other features of the projectile, such as excited or scattering states. 

There are a variety of state-of-the-art theoretical models to describe a two-body projectile breaking up in the 
Coulomb and nuclear fields of a target (fo r a review see 0). It should be noted that, although there are ways of 
dealing with three-body projectiles [13, 13 nere we wn l focus on two-body projectiles only. Regardless of the range 
of applicability, all these reaction models rely on an effective interaction between the two bodies that constitute the 
nucleus under study. Typically, the geometry of this interaction is fixed (a Woods-Saxon potential with tq ss 1.2 fm 
; 1 ■ and a ss 0.6 fm) and the depth is fitted to the binding energy of the projectile taking into account the correct 
angular momentum and possibly a spin-orbit force. It is clear that this potential is by no means unique and the 
uncertainty on the potential is expected to have an impact on the normalization of the breakup cross section. Within 
an effective-range theory, Typel and Baur have indeed shown that the electromagnetic strengths of halo nuclei depend 
on properties of the two-body projectile description These include bound-state properties, like the ANC, 

and parameters that characterize the final state interaction, like the scattering lengths. In Ref. it has also been 
found that the transfer cross section depend significantly on the single particle parameters. This dependence leads 
the authors to review the usual method for extracting spectroscopic factors from transfer reactions. 

Since the two-body system is loosely bound, one may think that the cross section should scale with the square of 
the ANC. Many questions remain open: is the reaction sufficiently peripheral for a dependence on ANC? To what 
extent is this true when couplings in the continuum are important? Or, how does the final state interaction affect 
the results? In this work we will look at Coulomb breakup from several angles: we study a proton halo 8 B and a 
neutron halo 11 Be, and look at the low energy regime (< 10 MeV/u) versus the higher energy regime (50-100 MeV/u). 
For this purpose, we use the Continuum Discretized Coupled Channel method |l8j or the time-dependent technique 
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|19|, according to applicability and feasibility. Other reaction models like DWBA or first order semiclassical theory 
|2fl EH ] are also used as qualitative tools to understand the results. In Sec. [HJ we present a brief description of the 
reaction models to be used. In Sec. IIIII the interactions describing realistic cases are specified as well as other details 
concerning the calculations. Sec. IIVI contains the results of the calculations of the Coulomb breakup of 8 B and 11 Be. 
In Sec. El tests on the sensitivity to the single particle parameters of the projectile are presented and discussed. 
Finally in Sec. IVII conclusions are drawn. 



II. THEORETICAL FRAMEWORK 



A. Projectile description 

Within the theoretical description of the breakup of a halo nucleus P impinging on a target T, the projectile is 
usually assumed to have a two-body structure: a pointlike and structureless fragment / (of mass to/ and charge Zje) 
loosely bound to a structureless core c (of mass m c and charge Z c e). The target is described as a structureless particle 
of mass tot and charge Ztg. The internal structure of the two-body projectile is described by the Hamiltonian Ho 

H = -^-A + V cf (r), (1) 

where r is the relative coordinate of the fragment to the core, and fi = m c mf/(m c + my) is the reduced mass of 
the c-f system. The potential V c f simulates the interaction between the core and the fragment. It is composed of a 
Coulomb term plus a nuclear term, which comprises a central part and a spin-orbit coupling term 

V cf (r) = V c (r, Rc) + V (r)+l- IV u (r), (2) 

where I is the orbital momentum of the c-f relative motion, and J is the spin of the fragment. The spin of the core 
is neglected and assumed to be nil. 

The Coulomb term Vc corresponds to the potential due to a uniformly charged sphere of radius Rq (the core) 
acting on a pointlike particle (the fragment) 




Vru. H,-)= { I /\ V r<Rc (3) 

r>R c . 

The central part of the nuclear potential Vb has a Woods-Saxon form factor 

V (r) = -Vif(r,R ,a), (4) 

where 

'r-Ro 



/(r, i? ,a) 



1 + exp 



a 

The spin-orbit coupling term has the usual Thomas form factor 



(5) 



Vu{r) =Vi S --^f(r,Ro, a). (6) 
r dr 

1 /3 

The radius is parameterized as: Ro = tqA c , where A c is the mass number of the core. The depths of the potential 
are adjusted to reproduce the bound states of the system and some of its resonances. 
In the Ij partial wave, the eigenstates of H are 

H (f> ljm (E,r) = Ecf> ljm (E,r), (7) 

where j is the total angular momentum resulting from the coupling of the orbital momentum I and the fragment spin 
J, and to is its projection. The negative-energy states correspond either to physical bound states or to Pauli-forbidden 
states simulating the presence of fragment-like particles in the core. They are normed to unity. The positive energy 
states describe the continuum of the two-body projectile, i.e. the scattering of the fragment by the core. Their radial 
part r~ l uij is normalized as 

uij(k,r) — >[cos5tj(k)Fi(k,r) +sin8ij(k)Gi(k,r)], (8) 



where k — y 2[iE/h 2 is the wave number of the relative motion of the core and the fragment, F; and G; are respectively 
the regular and irregular Coulomb functions |22| . and 81 j is the nuclear phase-shift. 
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B. CDCC 



The reaction of a loosely bound two-body projectile on a target can be approximated to a three-body scattering 
problem. Typically, target excitation is neglected and the three-body wave function is expanded in terms of the 
intrinsic motion of the core and the fragment c + / within the projectile uij(k, r), and the motion of the projectile 
relative to the target fhjj(k, R)- Here R is the vector connecting the center of mass of the projectile with the target, 
L is the corresponding orbital angular momentum, and J is the total angular momentum of the system. 

For modeling breakup it is important to have a good description of the projectile continuum. One of the most 
successful methods for elastic breakup is the Continuum Discretized Coupled Channel method (CDCC) [23| first 
developed for deuteron breakup but today widely applied to exotic nuclei (e.g. 0, 0> HE HE HE HI])- In the 
standard CDCC, the projectile continuum is included through a discretization into energy bins. For the bin wave 
function, we use the integral over the momentum interval, weighted by the function gij in the following way : 



„ / 2 f ki+1 

u ii,»( r ) = W^- J k gij{k)uij{k,r)dk . (9) 

For further details on the normalization iVy of the bins see Ref. [24| • 

The three-body Schrodinger equation can then be reduced to a coupled channel equation in R, 

f a j(R) = lL '~ L V a J :a ,(R)f a 'j(R), (10) 

where the subscript a represents all relevant quantum numbers and Ei is the average energy of continuum bin [ki, 
(or E < for the bound states), and £ is the total energy of the system. The coupling potential V^. a ,(R) consists 
of the sum of the core-target and fragment-target interactions averaged over the projectile states [29j . One can solve 
Eq. (|10[1 exactly: the solution \^ CDCC includes all possible rearrangements within the projectile continuum. One can 
also solve Eq. 1)1 iteratively: the one-step iteration of Eq. IjlOjl is equivalent to the standard DWBA, where the 
optical potential between the projectile and the target is replaced by the folding potential V J a , (R) j3(j ■ 



C. Time-dependent model 

The time-dependent description of the reaction 0, |^ HE HE HE HEl relies on the semiclassical approximation 
[2EEI in which the projectile-target relative motion is treated classically, while the internal motion of the projectile 
is described quantum mechanically. The target is assumed to follow a classical trajectory in the projectile rest frame. 
The projectile is therefore seen as evolving in a time-dependent potential which simulates its interaction with the 
target. The wave function 'J describing its internal structure is solution of the following time-dependent Schrodinger 
equation 

a |U(r, t) = [H + V(r, t)] *(r, t), (11) 

where Hq [see Eq. describes the internal structure of the projectile, and V is the time-dependent potential that 
simulates the P-T interaction. The latter reads 

V{r,t) = V cT [rMt)} + V fT [r f T(t)} - ^±^^ , (12) 

where the time-dependent coordinate R describes the classical trajectory followed by the target in the projectile rest 
frame. It is a hyperbola, which can be approximated fairly well by a straight line at high energies. The vectors r c x 
and TfT correspond respectively to the core-target and fragment-target coordinates. In Eq. (|12f> . V c t and VfT are 
local potentials which model the interaction between the target and the projectile constituents. They comprise a 
Coulomb term and a short-range optical potential, which simulates the nuclear interaction. 

Eq. Ijllj l is solved numerically considering the projectile initially in its ground state <pi j omo °f energy Eq. We use 
the algorithm described in Ref. [lflj . and consider linear trajectories. For each trajectory, characterized by impact 
parameter 6, we deduce the breakup probability Pbu by projecting the output wave function ^(r, t — ► +00) onto the 
positive eigenstates of Hq, which describe the projectile after dissociation. The breakup cross section Cbu is then 
obtained by summing this probability over all impact parameters. 
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Potential 


Vi 


Vis 


a 


ro 




(MeV) (MeV fm 2 ) 


(fm) 


(fm) 


Tl 


45.23 


19.59 


0.40 


1.25 


T2 


44.98 


19.59 


0.52 


1.25 


T3 


44.47 


19.59 


0.60 


1.25 


T4 


43.50 


19.59 


0.70 


1.25 


T5 


42.28 


19.59 


0.80 


1.25 



TABLE I: Parameters of the 7 Be-p potentials [see Eqs. ©-©]. Note that Ro used in (J^J-JSll is parameterized as roAl^ 3 , and 

1/3 

the Coulomb radius is always kept constant at Rc = 1-3A C fm. 



Potential 


Eq p 3/2 b 0p 3/2 

(MeV) (fm" 1 / 2 ) 


Eq p \/2 r 0p l/2 

(MeV) (MeV) 


Tl 


-0.137 0.6477 


2.8 2.2 


T2 


-0.137 0.7008 


2.3 1.6 


T3 


-0.137 0.7410 


1.9 1.2 


T4 


-0.137 0.7959 


1.6 0.9 


T5 


-0.137 0.8554 


1.3 0.7 



TABLE II: Levels obtained with the five potentials listed in TableQ] The energies and widths are expressed in MeV. The ANC 
of the 0p3/2 ground state b 0p3 / 2 is displayed as well. 



At sufficiently high energies, the solution of Eq. (|llfl can be approximated by first-order perturbation theory 
|2fl l2ll | . In that approximation the breakup process is assumed to occur in one step from the initial ground state to 
the continuum, and the couplings inside the continuum are neglected. The first-order approximation of the breakup 
cross section reads 



da 



(i) 



dE 



HE) 



1 2tt 
¥2j + 



T EE 

mo Ijm 



4>ij m {E,r) 



l V{r,t)dt 



db, 



(13) 



where lu — (E — E )/h. 

When the projectile-target interactions are purely Coulomb, the first-order approximation of the breakup probability 
can be easily expressed as a series of the contributions of each multipole of the time-dependent potential fll) . In 
practice, only the first terms of that series are needed. For straight-line trajectories, a semi-analytical expression of 
the dipole and quadrupole contributions to the first-order breakup probability can be obtained (see e.g. Ref. |37|). 



III. TWO-BODY INTERACTIONS 
A. Description of 8 B 



The 8 B nucleus is described by the usual two-body system: a proton loosely bound to a 7 Be core in its I ground 
state 0) ES EM US E3 • The internal structure of the core is neglected, and its spin is set to zero in the calculations. 
The 2 + ground state of 8 B is assumed to be a pure j?3/2 proton single-particle state. As explained in Sec. Ill Al the 
interaction between the core and the proton is simulated by a local potential [see Eqs. 

In order to study the sensitivity of the calculations to the description of the projectile, five 7 Be-p potentials are 
considered. The parameters of these potentials are listed in Tabled They are obtained by varying the diffuseness of a 
simplified version of the potential developed by Esbensen and Bertsch (potential T2). While the spin-orbit term 
is kept unchanged for all potentials, the depth of the central part Vi is adjusted to reproduce the 137 keV binding 
energy of the bound state [4l|. Since the spin of the core is neglected, the + , 1 + , and 3 + states resulting from the 
coupling of this | spin and the | angular momentum of the halo proton are degenerate with the 2 + ground state. 
Therefore, the p-wave 1 + resonance at E = 0.63 MeV |4l| is not reproduced by this parameterization. We use the 
same potential for all partial waves. The levels obtained with the potentials of Table [I] are listed in Table [H] The 
ANC of the ground state &op3/2 is given as well. Note that besides the 0p3/2 ground state, the potentials exhibit 
also a pl/2 resonance, which does not correspond to any known physical state. This unphysical state has not been 
adjusted, and its energy varies from one potential to the other. 
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Potential V; eV en V( od d Vis a r 

(MeV) (MeV) (MeV fm 2 ) (fm) (fm) 



TABLE III: Parameters of the 



VI 


62.52 


39.74 


21.0 


0.6 


1.2 


V2 


66.325 


38.37 


12.44 


0.5 


1.2 


V3 


58.905 


40.025 


27.68 


0.7 


1.2 


V4 


71.28 


49.015 


29.95 


0.6 


1 1 
1 . 1 


V5 


55.25 


32.515 


12.86 


0.6 


1.3 


V6 


59.05 


59.05 





0.62 


1.236 


n potentials [see Eqs. Note that R 


o used 


Potential 






Eopl/2 


-E'0d5/2 


rod5/2 




(MeV) (fm" 1/2 ) 


(MeV) 


(MeV) (MeV) 


VI 


-0.504 


0.83 


-0.184 


1.274 


0.162 


V2 


-0.504 


0.80 


-0.184 


1.274 


0.131 


V3 


-0.504 


0.87 


-0.184 


1.274 


0.200 


V4 


-0.504 


0.82 


-0.184 


1.274 


0.145 


V5 


-0.504 


0.85 


-0.184 


1.274 


0.181 


V6 


-0.504 


0.85 


-12.5 


2.6 


0.9 



used in ©-© is parameterized as tqA, 



1/3 



TABLE IV: Levels obtained with the six potentials listed in Table lllll Potentials VI to V5 were adjusted to the first three 
levels of 11 Be, while V6 reproduces only the | + ground state. The ANC of the lsl/2 state is displayed as well. 



B. Description of n Be 

As in previous studies 0, 0, H^, |3 HE E3) E3> 01 > 11 Be is described as a neutron loosely bound to a 10 Be core. 
The 10 Be core is assumed to be in its + ground state. The depths of the V c f potential J5J simulating the interaction 
between 10 Be and the neutron are adjusted to reproduce the low-lying levels in the 11 Be spectrum 45]. The well 
known shell inversion observed between the two bound states of 11 Be is reproduced using a parity dependent depth of 
the central term of the nuclear potential Vi. The i + ground state is modeled by a lsl/2 state, the \ excited state 

by a Opl/2 state, and the | + resonance is adjusted in the d.5/2 partial wave. 

In this case, we make use of six sets of parameters to study the sensitivity of our calculations to the potential choice. 
They are listed in Table 1TTT1 The first potential (VI) was developed for a time-dependent analysis of the breakup 
of 11 Be on 12 C 0|. The next four (V2 to V5) have been derived from VI by varying either the diffuseness or the 
radius of the Woods-Saxon form factor The values were chosen to encompass those used in most other breakup 
calculations [3 ll 1321 l36l l42j| . In addition, we also use a sixth potential (V6) developed by Fukuda et al. for analyzing 
their data [Jit reproduces only the ground state energy of 11 Be and does not include a spin-orbit coupling term. 

The physical energy levels obtained with these six potentials are listed in Table IIVI The asymptotic normalization 
coefficient b lsl / 2 of the ground state is also given. 



C. Projectile-target interactions 

The reactions we consider are Coulomb dominated. However, the nuclear interactions between the target and the 
projectile composites are not negligible |l9l l24|. As mentioned earlier, they are simulated by short ranged optical 
potentials. The general form factor of the potential simulating the interaction between the constituent x and the 
target T is 

V xT (r) = V c {r, Rc) - Vf(r, R R , a R ) - iWf(r, Ri, ai ) - iW Dai -^-f{r, Ri, ai ), (14) 

ar 

where Vc is the point-sphere Coulomb potential ©, and / is the Woods-Saxon form factor (JSJ. The values of the 
parameters we consider for these interactions are listed in Table E] They correspond to potentials found in the 
literature for similar systems. 

As in Ref. l24l . we simulate the interaction between 7 Be and 58 Ni at low energy using the potential developed by 
Moroz et al. [43. This potential reproduces elastic scattering data of 7 Li impinging on 58 Ni at 14.2 MeV. 

To model the interaction between 10 Be and 208 Pb at 69 MeV per nucleon, as in Ref. we follow Typel and 
Shyam an d adapt a potential developed by Bonin et al. , which reproduces the elastic scattering cross section 
of a on lead at 699 MeV. 
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(MeV) 




Cfml (MeV) (Me VI 

I 1111 / 1 1V1C V / \ 1V1\_^ V / 




(fm) 




'Be 


as Ni 


3 


100 


4.064 


0.65 30.6 


4.347 


0.80 


5.032 


10 Be 


208p b 


69 


70.0 


7.43 


1.0 58.9 


7.19 


1.04 


5.92 


P 


58 Ni 


3 


54.512 


4.529 


0.75 11.836 


4.877 


0.58 


5.032 


n 


208p b 


69 


29.46 


6.93 


0.75 13.4 


7.47 


0.58 





TABLE V: Parameters of the optical potentials 1141 which simulate the interaction between the projectile constituents and the 
targets. 

The interactions of the valence nucleon with the targets are approximated by the parameterization of Becchetti and 
Greenlees |Z^, as m previous works [ilil24j| . 

IV. BREAKUP CALCULATIONS 
A. The 8 B case 

The breakup of 8 B on 58 Ni has been measured at 25.75 MeV at Notre-Dame [53,|^. This reaction is analyzed in 
Ref. |24j within the CDCC framework (see Sec. Ill B|) . The calculation is in good agreement with the data using the 
two-body description of 8 B given in Sec. EEI (potential T2 of Table 0. 

In the present paper, we study the sensitivity of this analysis to the 7 Be-p potential choice. CDCC calculations 
identical to that of Ref. are performed with the code FRESCO [5^ considering the five potentials T1-T5 listed 
in Table [I] For the model space, we use L up to 1000, a maximum radius for the distorted waves of i? m ax = 500 fm 
and coupling potentials truncated at i?bin = 60 fm. The discretization of the continuum of the projectile is made up 
to E max = 4 MeV and all dipole and quadrupole transitions for s, p and d waves are included. As in Ref. the 
nuclear interactions between the projectile constituents and the nickel target are simulated by the optical potentials 
given in Table Ivl 

First, let us consider simpler cases of the continuum wave functions whilst changing the ground state interaction 
according to models T1-T5. If only pure Coulomb waves are taken in the 7 Be-p continuum, the scaling with the 
square of the ANC of the initial bound state b 0p3 / 2 is perfect (i.e. with less than 1 % difference). This same result 
is obtained if the nuclear scattering 7 Be-p potential is kept fixed to Tl for all models. We have checked these results 
using other reaction models. The DWBA approximation leads to the same conclusion, as well as the first-order 
approximation when the breakup of 8 B on 58 Ni is computed at higher energy (50 MeV/nucleon) considering a purely 
Coulomb interaction between the projectile and the target. The direct proportionality of the cross section to the 
square of the ANC of the initial bound state indicates that only the asymptotic part of the wave function has an 
influence upon the cross section, consistent with the small binding energy of the projectile. 

When the scattering potential is chosen to correspond to the ground-state interaction, the breakup cross section 
is only approximately proportional to b^ p3 ^ 2 . This is illustrated in Fig.^ where the breakup cross sections obtained 

with potentials T1-T5 are plotted as a function of the relative energy between the 7 Be core and the proton after 
breakup. To emphasize the sensitivity to the potential choice, those values are divided by the square of the ANC. 
The different curves exhibit similar behaviors. They all increase sharply to reach a maximum around 0.5 MeV, and 
display a bump in the slow decrease that follows the maximum. However, these behaviors are not identical. First, 
even though they have been divided by bQ p3 , 2 , their magnitudes differ by approximately 4 % in the vicinity of the 
maximum. Second, the location of the bump changes with the potential choice. Since these features do not appear 
when the same scattering potential is used, they are due to differences in the distorted waves. 

To understand these features, we analyze the contributions of the different partial waves to the breakup cross 
section. The sensitivity of these contributions to the potential choice vary from one partial wave to the other. This is 
illustrated in Fig. [2 which displays the p3/2 (dominant) and the pl/2 contributions to the cross section. They both 
have been divided by &Q p3 / 2 - O n the one hand, the p3/2 contributions — as well as the s and d ones — exhibit very 
similar shapes, but their magnitudes vary by 4 % — this variation goes up to 8 % for the s component. These results 
explain the non exact proportionality of the total breakup cross section to the square of the single particle ANC. On 
the other hand, each pl/2 contribution exhibits a maximum whose location varies with the potential choice. These 
maxima are responsible for the presence of the bumps in the total cross sections above 1 MeV. They are related to the 
pl/2 resonance obtained with potentials T1-T5 (see Table ITTj l. This effect has already been observed in the breakup of 
11 Be on 12 C, where the presence of the | + resonance in the 11 Be spectrum induces a narrow peak in the cross section 
|fj.l44|. It confirms that the description of the continuum has a significant influence on the breakup calculation. 
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FIG. 1: Cross section for the breakup of 8 B on 58 Ni at 25.75 MeV divided by the square of the single particle ANC. The cross 
section is given as a function of the 7 Be-p relative energy E after breakup. Calculations are performed with the different 7 Be-p 
potentials listed in Table [I] 




0.5 1 1.5 2 2.5 3 

E (MeV) 



FIG. 2: Contributions of p3/2 and pl/2 partial waves to the breakup cross section of 8 B on 58 Ni. They are plotted as a function 
of E, and are divided by the square of the single particle ANC. Calculations are performed with the different 7 Be-p potentials 
listed in Table Q 



The same effects are observed in DWBA calculations, and, at higher energy (50 MeV/nucleon), within the first-order 
perturbation theory. However, since the contribution of the E2 component decreases at high energy, the sensitivity 
to the p-wave continuum is less noticeable in the latter case. 

These results indicate that the breakup calculations are sensitive, not only to the tail of the ground-state wave 
function, but also to the way the continuum of the projectile is described. This influence of the continuum description 
upon breakup calculations has already been mentioned in Ref. |Tll |. In that reference, Trache et al. observed that the 
momentum distributions computed using distorted waves differed from those obtained previously with plane waves 
[Tof. Since the reaction is very peripheral, the sensitivity to the scattering potential is most likely due to the subsequent 
variations in the phase shifts. The detailed analysis of these effects is addressed in Sec. El 
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FIG. 3: Breakup cross section of n Be on 208 Pb at 69 MeV/nucleon as a function of the energy E. The value is divided by the 
square of the single particle ANC. Calculations are performed with potentials V1-V6 listed in Table IV1 



A more severe case of the breakdown of the ANC scaling is found in the Coulomb breakup of 11 Be on lead. This 
reaction has recently been measured at RIKEN at 69 MeV/nucleon 0. From the analysis of their experiment, Fukuda 
et al. extracted a spectroscopic factor of about 0.7 for the lo Be(0 + ) ®sl/2) configuration of the ground state of n Be. 

With the aim of analyzing the sensitivity of this figure to the 10 Be-n potential, we perform time-dependent calcu- 
lations of the reactions using the six potentials V1-V6 described in Sec. IIII 51 As in Ref. |19|. the calculations are 
done within the time-dependent framework described in Sec. Ill CI using the optical potentials listed in Table for 
simulating the nuclear interactions between the projectile and the target. 

As for 8 B, if simplifications are performed in the continuum (namely, switching off the nuclear Be-n interaction, or 
taking the same interaction for all scattering partial waves as the ground state) , the resulting breakup cross section is 
directly proportional to the square of the ground state ANC b\ sX , 2 . When the realistic 11 Be interactions of table II are 
used, the situation changes: the cross section is no longer proportional to b\ sl i 2 - This is illustrated in Fig. where 
the breakup cross sections divided by b\ sl , 2 are plotted as a function of the relative energy between the 10 Be core and 
the neutron after breakup. All curves exhibit the same shape — if one excepts the small bump around 1.3 MeV due to 
the db/2 resonance, which is not reproduced by the potential V6. However, their magnitude varies significantly from 
one potential to the other. In particular, V6 leads to a cross section larger by 40 % than V1-V5. Moreover, even 
though they have been adjusted on the same energy levels, these five potentials lead to variations in the breakup cross 
sections as large as 20 %. Since this discrepancy is observed only when the scattering potential differs, we conclude 
that the way the continuum is described has a significant influence on the breakup calculation. The same result is 
obtained within the first-order approximation. With such a variation of the calculation with the potential choice, the 
validity of the spectroscopic factor extracted by Fukuda et al. from their breakup measurement |!j is questionable. 

The effect observed in Fig.[3]is very different from one partial wave to the other. This is illustrated in Fig.QJ where 
the dominant p3/2 and pl/2 contributions to the total cross section are depicted. Like the results shown in Fig. [31 they 
have been divided by b\ sl j 2 to cancel the dependence on the ANC. On the one hand, the difference between V1-V5 
is located mainly in the p3/2 contribution. These potentials lead indeed to very similar j?l/2 contributions. On the 
other hand, even if the j?3/2 contribution of potential V6 is still higher than that of the others, the major difference 
between that potential and the other five lies mainly in the pl/2 contribution. This partial wave analysis confirms 
that the breakup cross section depends not only on the ground-state wave function, but also on the description of the 
partial waves describing the continuum. 



B. The "Be case 
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FIG. 4: Contributions of the main partial waves p3/2 and pl/2 to the breakup cross section of n Be on 208 Pb. The values are 
divided by the square of the single particle. 

V. ANALYSIS OF THE SENSITIVITY TO THE POTENTIAL CHOICE 

The features presented in the previous section are all reproduced using fully coupled reaction models (like CDCC or 
time-dependent approaches) and also one-step approximations (like DWBA or first-order perturbation theory). This 
indicates that those dependences appear mainly through first-order transitions from the initial bound state to the 
continuum. Therefore, higher-order effects will be neglected in the qualitative analysis presented in this section. The 
effects observed in section llVl can be explained through the study of the wave functions of the initial bound states 
and those of the partial waves describing the continuum [see Eq. Since for both nuclei the breakup reaction is 

Coulomb dominated, we neglect the nuclear contribution in the following qualitative analysis and consider only the 
Coulomb potential, which we expand into multipoles. 

The contributions of the p waves to the breakup of 8 B on 58 Ni displayed in Fig.Elcorrespond mainly to E2 transitions 
from the ground state. These transitions are large because the reaction occurs at low energy. Fig. pictures the 
radial wave functions of the ground state (a) and the continuum states, p3/2 (b) and pl/2 (c), computed at an energy 
of 1 MeV. 

Once divided by their ANC, the wave functions of the 8 B ground state obtained with the different 7 Be-p potentials 
are identical above 5 fm [Fig.OJa)]. The only difference between them lies below that radius. Due to the low binding 
energy of the system, these wave functions exhibit a very long range. This long range, and the presence of the r 2 
factor of the quadrupole term of the Coulomb potential in the transition matrix element ensures the periphcrality 
of the breakup reaction. Consequently, the breakup cross sections are not very sensitive to the differences at small 
radius. This explains why, when the same description of the continuum is chosen in all calculations, the cross section 
is exactly proportional to &op3/2' ^ the continuum is described using potentials T1-T5, the corresponding wave 
functions differ from one potential to the other. Since the reaction is very peripheral, only the differences in the 
distorted waves for r > 5 fm (i.e. mainly phase shifts) lead to variations in the breakup cross sections. This difference 
is illustrated in part (b) of Fig.[S]for the dominant p3/2 partial wave. Potentials T1-T5 indeed lead to variations in 
the phase shift of the wave function. These remain small though due to the fact that all potentials have been adjusted 
to reproduce the low-lying ground state in that partial wave. However, they are not negligible, and explain the loss 
of exact proportionality of the p3/2 contribution to the breakup cross section to the square of the ANC, illustrated 
in Fig. |21 A similar analysis of the s and d continuum wave functions explains the lack of proportionality to &o P 3/2 °^ 
these components as well. 

The differences between the pl/2 partial waves are much more significant [Fig. Etc)]. This is due to the presence in 
that partial wave of a low- lying resonance whose energy varies from one potential to the other (see Table |nj. These 
large variations account for the significant distortions shown in Fig. [2] 

In the case of the Coulomb breakup of 11 Be at 69 MeV/nucleon the same reasoning can be done in order to explain 
the significant variations observed in Fig. |3J Occurring at higher energy, the transition to the continuum is dominated 
by the El term of the Coulomb interaction |37J. Therefore, starting from an initial s state, the transfer to the 
continuum occurs mainly through p waves. 
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FIG. 5: Wave functions of 8 B obtained with the different 7 Be-p potentials T1-T5: (a) the initial ground state divided by &op3/2, 
(b) and (c) the p3/2 and pl/2, respectively, continuum wave functions computed at E — 1 MeV. 

Fig. depicts the radial wave functions of the initial lsl/2 ground state (a), and of the p3/2 (b) and pl/2 (c) 
scattering states at E = 1 MeV. They are calculated for the potentials V1-V6 of Table ITTT1 As in Fig. Eta), the 
ground state wave function is divided by the ANC. As in the 8 B case, due to the long range of the wave function, 
and to the r factor appearing in the matrix element, the first order breakup cross section is not sensitive to the 
variations at small radii. The cross section should therefore be proportional to b\ sl / 2 - However, as for 8 B, the change 
in the potential leads also to variations in the scattering wave functions. In the present case, though, the changes 
in the dominant component are much more significant than in the previous example, explaining the large differences 
observed in the breakup cross section in Sec. HVBl We indeed observe significant variations in the p3/2 phase shift 
[Fig. Efb)]. These variations are responsible for the large differences in the p3/2 contribution to the breakup cross 
section observed in Fig. 0] They overcome by far the differences due to the ANC, which means that the 11 Be breakup 
reaction of Ref. seems more sensitive to the way its continuum is described than to the details of its initial single 
particle state. 

The differences observed in the pl/2 contribution to the cross section (Fig. 0} can easily be understood by the 
variations in the phase shift shown here. The very similar pl/2 phase shifts obtained with potentials V1-V5 justify 
that these five potentials lead to approximately the same pl/2 contributions to (da/d^l)/b^ sl ^ 2 . The similarity in 
the pl/2 phase shifts is due to the fact that all those potentials have been adjusted in order to reproduce the very 
loosely bound \ excited state of 11 Be in that partial wave. Potential V6, that does not reproduce this excited state, 
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FIG. 6: Wave functions of 11 Be obtained with the different 10 Be-n potentials V1-V6: (a) the initial ground state divided by 
bisi/2, (b) an d (c) the p3/2 and pl/2, respectively, continuum wave functions computed at E — 1 MeV. 

presents a very different phase shift of the pl/2 component, which explains its larger pl/2 contribution to the cross 
section. 

This first-order analysis of the breakup reactions explains qualitatively the sensitivity of our calculations to the 
core-fragment potential. It shows that the wave functions describing both the ground state and the continuum have a 
significant influence on the cross section. Since the reaction is mostly peripheral, this influence occurs mainly through 
the asymptotic characteristics of the wave functions, i.e. the ANC of the initial bound state, and the phase-shift in 
the continuum. 

These results confirm the study of the electromagnetic strengths in one-nuclcon halo nuclei performed by Typel and 
Baur in Ref. [16|. In that paper, the authors analyze the sensitivity of i?(EA) transitions to the two-body description 
of the loosely-bound nuclei. They show that the energy distributions depend not only on the description of the bound 
state, but also on that of the continuum. Assuming the actual wave functions can be replaced by their asymptotic 
behaviours, they obtain simple analytical expressions for the energy distributions. These expressions show clearly that 
the transition strengths can be characterized by only a few parameters describing the nucleus: the binding energy of 
the system, the ANC of the bound state wave function, and the scattering length describing the final state interaction. 
The present analysis shows that these conclusions remain qualitatively valid in breakup reactions where higher-order 
effects and nuclear interactions are significant [2J, |3jj • 

It should be noted that the shape of the transition strengths obtained by Typel and Baur depends significantly 
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on the scattering length (see Figs. 9 and 10 of Ref. |Tq). If this were so, it should enable the extraction of the 
ANC and the scattering length from the same data. However, our results show hardly any distortion of the breakup 
cross sections with the potential choice. The curves differ only in amplitude, not in shape, even though the various 
potentials lead to significantly different phase shifts. The consequence is then that the final state interaction needs 
to be constrained by other observables in addition to the breakup energy distribution. This discrepancy between our 
results and those of Ref. is due to further approximations performed in Ref. required for obtaining analytic 
solutions. We find that the analytic expressions of Ref. are only useful qualitatively. The shape and peak of the 
energy distributions are not correctly produced. 

VI. CONCLUSION 

In this paper, the breakup calculations of two loosely bound projectiles, 8 B and n Be, have been performed within 
CDCC and time-dependent frameworks at low and high energy, respectively. Different potentials were used to simulate 
the interaction between both constituents of the projectile. It has been shown that the potential geometry has a 
significant influence upon the breakup cross section. This dependence occurs through the asymptotic normalization 
coefficient of the initial bound state, as expected from the peripheral nature of the reaction. However, the way excited 
states and the continuum are modeled also plays a role. 

These results have been confirmed using DWBA and/or first-order perturbation theory. These approximations 
enabled us to interpret the variations observed in the more elaborated techniques with one-step theory. Our analysis 
shows with transparency that different potential geometries lead to differences in the phase shifts that can affect 
significantly the total cross section. This qualitative analysis is in good agreement with Typel and Baur's analytical 
study of the electromagnetic strengths in halo nuclei . 

One of our test cases consisted of the breakup of 8 B on 58 Ni at 25.75 MeV. Different potentials lead to minor 
variations in the dominant p3/2 phase shift, since all potentials have been adjusted to reproduce the 8 B ground state 
in that partial wave. These variations induce only 4 % differences in the breakup cross section, which remains therefore 
approximately proportional to the square of the ANC. On the contrary, our potentials give large differences in the 
pl/2 phase shift, leading to significant distortion in the energy distribution for this partial wave. As this partial wave 
is less important in this reaction, the phase shift differences have a smaller effect in the normalization of the total 
cross section. 

In the breakup of n Be on 208 Pb at 69 MeV/nucleon, no fitting was imposed on the dominant p3/2 partial wave, 
thus the various potentials lead to significant differences in the corresponding phase shifts. These differences are so 
large that they overcome those in the single particle ANC. The cross section is therefore not proportional to the 
square of the ANC. Five out of the six potentials we have used for this analysis have been adjusted to reproduce 
the \ excited state of n Be in the pl/2 partial wave. This induces very similar phase shifts in that partial wave. 
That contribution to the cross section is therefore nearly proportional to the square of the ANC. A sixth potential, 
that does not reproduce the excited state, has a very different pl/2 phase shift, and leads to a much higher pl/2 
contribution. 

These results show that Coulomb breakup reactions probe not only the ground state of the projectile, as it is usually 
assumed, but also the continuum description. We have indeed seen that differences in the phase shift of the dominant 
partial wave can influence significantly breakup calculations. Therefore, one should be particularly cautious when 
extracting ground-state spectroscopic information or ANCs from breakup measurements. As shown in the calculation 
of the Coulomb breakup of 11 Be, different descriptions of the continuum can lead to variations in the cross section 
up to 40 %, generating a significant inaccuracy in the spectroscopic factor. Moreover, even though the reaction is 
very peripheral, the extraction of the single-particle ANC from breakup measurements can be very tricky due to the 
strong dependence of the calculations on the scattering waves. The sensitivity of the extracted values to the potential 
geometry should therefore be evaluated, and taken into account in the analysis. Whenever possible, the potentials 
used to simulate two-body projectiles should be adjusted on other experimental data (excited states and phase shifts). 
This underlines the need for additional data to constrain the scattering properties of the core-fragment potentials. 
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